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Abstract 

In this paper the Boltzmann transport equation for thermal processes induced by ultra-short 
laser pulses is formulated and solved. For thermal process where the duration of the laser 
pulse (femtoseconds - attoseconds) is shorter than the relaxation time, the solution of 
Boltzmann equation differs from the solution of the Fourier equation. In this paper the 
formula for density current of the heat carrier is obtained. It is shown that for thin one- 
dimensional structures the current strongly depends on the scattering mechanism of the heat 
carriers. 
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1. Introduction 

Recently it has been shown that after optical excitation by femtosecond pulse establishment of 
an electron temperature by e-e scattering takes place on a few hundred femtosecond time 
scale in both bulk and nanostructured noble materials [1-4]. In noble metal clusters the 
electron thermalization time (relaxation time) is of the order of 200 fs [3,4]. This relaxation 
time is much larger than the duration of the now available femtosecond optical pulses offering 
the unique possibility of analyzing the properties of a thermal quasi-free electron gas [5]. In 
paper [5] using a two color femtosecond pump-probe laser technique the ultrafast energy 
exchanges of a nonequilibrium electron gas was investigated. When the duration of the laser 
pulse, 25 fs, in paper [5], is shorter than the relaxation time the parabolic Fourier equation 
cannot be used [6,7]. Instead, the new hyperbolic quantum heat transport equation is the valid 
equation [8]. The quantum heat transport equation is the wave damped equation for heat 
phenomena on the femtosecond scale. 

Wave is an organized propagating imbalance [10]. Some phenomena seem to be clearly 
diffusive, with no wave-like implications, heat for instance. That was consistent with 
experiments at the late century, but not any longer. As far back as the 1960s ballistic (wave- 
like) heat pulses were observed at low temperatures [10]. The idea was that heat is just the 
manifestation of microscopic motion. Computing the classical resonant frequencies of atoms 

13 

or molecules in a lattice gives numbers of the order of 10 Hz, that is in the infrared, so when 
molecules jiggle they give off heat. These lattice vibrations are called phonons. Phonons have 
both wave-like and particle-like aspects. Lattice vibrations are responsible for the transport of 
heat, and we know that is a diffusive phenomenon, described by the Fourier equation. 
However, if the lattice is cooled to near absolute zero, the mean free scattering of the phonons 
becomes comparable to the macroscopic size of the sample. When this happens, lattice 
vibrations no longer behave diffusively but are actually wave-like or thermal wave. By 
controlling the temperature of a sample, one can control the extent to which heat is ballistic 
(thermal wave) or diffusive. In essence if a heat pulse is launched into sample (by the laser 
pulse interaction) and if the phonons can get across the sample without scattering, they will 
propagates as thermal waves. 

The extent to which the motion of quasiparticles (phonons) or particles is ballistic, is 
described by the value of the relaxation time, x . For ballistic (wave-like) motion, x — > °° . 
The equation which is generalization of the Fourier equation (in which x — > °° ) is the 
Heaviside equation [6] for thermal processes: 

d 2 T dT 

For very short relaxation time, T — > we obtain from equation (1) the Fourier equation 



+ — = DV 2 T (1) 
dt 

from equation (1) the Fourier equation 
= DV 2 T (2) 



dt 

and for T — > °° we obtain from formula (1), the ballistic = thermal wave motion: 



1 d 2 T 

2 dt 2 



= v 2 r. (3) 



v 

In the set of papers [6,7,8] the quantum generalization of the Heaviside equation was obtained 
and solved: 

1 d 2 T | m dT _ d 2 T 
v 2 dt 2 ri dt dt 2 

where v = ccc and is the fine structure constant, c is the vacuum light velocity. In formula (4) 
m is the heaton mass [3]. Heaton energy is equal 
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E h =ma 2 c 2 . (5) 

In papers [6] Heaviside equation was obtained for the fermionic gases (electrons, nucleons, 
quarks). In this paper the Heaviside equation will be obtained for particles with mass m, 
where m is the mass of the fermion or boson. Moreover beside the elastic scattering of the 
particles, the creation and absorption of the heat carriers will be discussed. The new form of 
the discrete Heaviside equation will obtained as the result of the discretization of the one- 
dimensional Boltzmann equation. The solution of the discrete Boltzmann equation will be 
obtained for Cauchy boundary conditions, initialed by ultra-short laser pulses, i.e. for 
At < T, the relaxation time. 

2. The model 

Let us consider the one-dimensional rod (strand) which can transport "particles" - heat 
carriers. These particles, however may move only to the right or to the left on the rod. Moving 
particles may collide with the fixed scatter centra, barriers, dislocations) the probabilities of 
such collisions and their expected results being specified. All particles will be of the same 
kind, with the same energy and other physical specifications distinguishable only by their 
direction. 
Let us define: 

u(z,t) = expected density of particles at z and at time t moving to the right, 
v(z,t) = expected density of particles at z and at time t moving to the left. 
Furthermore, let 

8(z) = probability of collision occurring between a fixed scattering centrum and a particle 
moving between z and z + A. 

Suppose that a collision might result in the disappearance of the moving particle without new 
particle appearing. Such a phenomenon is called absorption. Or the moving particle may be 
reversed in direction or back-scattered. We shall agreeing that in each collision at z an 
expected total of F(z) particles arises moving in the direction of the original particle, B(z) 
arise going in the opposite direction. 

The expected total number of right-moving particles in z l < z < z 2 at time t is 

^u(z,t)dz, (6) 
while the total number of particles passing z to the right in the time interval t x <t<t 1 is 



vju(z,t)dt, (7) 



wj 

h 

where w is the particles speed. 

Consider the particle moving to the right and passing z + A in the time interval 

t.+y <t<t 2 +y 

1 /w 2 /w 



t 2 +t\l w 1 2 / A\ 

w ^u(z + A,t')dt' = w^u \z + A,t'+— it'. (8) 



w 

t, + A / w 1, V / 



These can arise from particles which passed z in the time interval t x <t<t 2 and came through 
( z, z + A ) without collision 

h 

w\(l-A8(z,t'))u(z,t')dt' (9) 
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plus contributions from collisions in the interval (z,z + A). The right-moving particles 
interacting in ( z, z + A ) produce in the time t\ to ? 2 , 



"2 

w§A8(z,t')F(z,t')u(z,t')df 

h 

particles to the right, while the left moving ones give: 



JjA8(z,t')B(zJ')r(zJ')ch' 



(10) 



(11) 



Thus 



w 



z + A,t'+— dt'= w\u(z,t')dt'+wA\d(z,t')(F(z,t')-l)u(z,t')dt' 

W J I i 

h 

+wA^d(z,f)B(z,f)v(z,f)dt'. 



(12) 



Now, we can write: 



to get 



( 



z + A,t'+- 



w 



u(z,t') + 



du 



1 du 



(z,n+~(z,n 

az w at 



(13) 



rf au< 13 



(z,f) +~(z,t')dt'= \5(z,t')((F(z,t')-lMz,t') + B(z,t'Mz,t'))dt'. 
dz ) wot 



On letting A — > and differentiating with respect to h we find 



In a like manner 



^+-^ = d(z,t)(F(z,t)-lMz,t) + 5(z,t)B(z,t)v(z,t). 
az w at 



-^ + —^ = 8(z,t)B(z,t)u(z,t) + 8(z,t)(F(z,t)-l)v(z,t). 
az w at 



(14) 
(15) 

(16) 



The system of partial differential equations of hyperbolic type (15,16) is the Boltzmann 
equation for one dimensional transport phenomena [11]. 
Let us define the total density for heat carriers, p(z,t) 

p(z,t) = u(z,t) + v(z,t) (17) 

and density of heat current 

j(z,t) = w(u(z,t)-v(z,t)). (18) 
Considering equations (15-18) one obtains 
dp , 1 dj 



dz + w 2 dt 



= 8(zj)u(zj)(F(z,t)-B(z,t)-l) + b(z,tMz,t)(B(z,t)-F(z,t) + l). 



Equation (19) can be written as 



dp { 1 dj = S(z,t)(F(z,t)-B(z,t)-l)j 
dz w 2 dt 



w 



or 



J = 



w 



dp 



+ - 



1 



S(z,t)(F(z,t)-B(z,t)-l) dz wb(z,t)(F(z,t)-B(z,t)-l) dt 



(19) 
(20) 

(21) 
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Denoting, D, diffusion coefficient 

w 



D = 



S(z,t)(F(z,t)-B(z,t)-l) 



and x , relaxation time 



equation (21) takes the form 



x = 1 (22) 

w8(z,t)(l-F(z,t)-B(z,t)) 

j = ~D d l-4- (23) 
dz dt 

Equation (23) is the Cattaneo's type equation and is the generalization of the Fourier 
equation. Now in a like manner we obtain from equation (15-18) 

~+~=S(z,t)u(z, t)(F(z, t)-l + B(z, t)) 

w oz w at (^4) 

+8(z,t)v(z,t)(B(z,t) + F(z,t)-l)) 

or 

|U|e=a (25) 

dz dt 

Equation (25) describes the conservation of energy in the transport processes. 

Considering equations (23) and (25) for the constant D and x the hyperbolic Heaviside 

equation is obtained: 

dt 2 dt dz 

In the case of the heaton gas with temperature T(z,t) equation (26) has the form [6,7,8] 

d 2 T dT _ d 2 T 

% dt 2 + dt ~ dz 2 ' 
where x is the relaxation time for the thermal processes. 

3. The solution of the Boltzmann equation for the stationary transport 
phenomena in one dimensional strand 

In the stationary state transport phenomena dF(z,t)/ dt = dB(z,t)dt = and dd(z,t)/ dt = 0. 
In that case we denote F(z,t) = F{z) = B(z,t) = B(z) = k(z) and equation (10) and (11) can 
be written as 

^=8(z)(k-l)u(z) + 8(z)kv(z), 

dz (27) 

- — = 8(z)k(z)u(z) + 5(z)(k(z) - l)v(z) 
dz 



i^r+^ = D-i. (26) 



with diffusion coefficient 



and relaxation time 



w 

D = — (28) 
5(z) 



x(z) = -= 1 • (29) 

wd(z)(l-2k(z)) 



The system of equations (27) can be written as 
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,2 — (Sfc) . 



' 8 »(2t-l) + *0-« + - 8<t - 1), ' (8ty 



8/: dz 



du 



= 8(k-l)u + 8kv(z). 



Equation (30) after differention has the form 

d 2 u ,. .du , . , . _ 
— + f(z) — + g(z)u(z) = 0, 



where 



dz l 



f(z) = 



dz 



(30) 
(31) 

(32) 



1 



8 dk dd 
+ — 

k dz dz 



g(z)=d 2 (z)(2k-l)- 
For the constant absorption rate we put 

k(z) - k - constant ^ 
In that case 



8 dk 



(33) 



k dz 
1 



f(z) = 



l_dS 
8 dz 



(34) 



g(z)=5 2 (z)(zk-l). 

With functions /(z) and g(z) the general solution of the equation (30) has the form 



u(z) = C x e ) ^ z +C 2 e 



(35) 

In the subsequent we will consider the solution of the equation (32) with f(z) and g(z) 
described by (34) for Cauchy condition: 

u(0) = q, v(a) = 0. (36) 
Boundary condition (36) describes the generation of the heat carriers (by illuminating the left 
end of the strand with laser pulses) with velocity q heat carrier per second. 
The solution has the form: 



u(z) 



2qe 



[/(0)-/(a)] 



1 + J g e 2[/(0)-/( fl )] 

k-\ 



(1-2*) 



cosh[/Xx)-/(a)] 



u(z) = 



(l-2k) 2 -(k-l) 

2qe (m - f(a)) 



(l-2k) 2 -(k-l) 
— sinh|/(x)-/(a)l 



(37) 



l+(3e 



2l/(0)-/( fl )J 



k 



where 
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(38) 



f(z) = {\-2ky\5dz, 

f{Q)={\-2ky[\5dz\, 

f(a) = (1-2*)* [\ddz\, 
i 

(\-2ky+(k-i) 
{i-iky ~{k-\) 

Considering formulae (17), (18) and (37) we obtain for the density, p( z ) and current 
density j(z). 

2qwe [f(Q) - f(a)] 



Kz) = 



2[/(0)-/(a)J 



(1-2*) 2 
i 

(1-2*) 2 -(*-l) 



-cosh[/-(z) -/(«)]- 



1-2* 



-sinh|/(z) -/(a)] 



(1-2*) 2 -(*-l) 



and 



2qe 



[/(0)-/(a)] 



1+jSe 



2l/(0)-/( fl )J 



(l-2*) : 



(1-2*) 2 -(*-l) 



-cosh[/-(z) -/(«)]- 



(1-2*) 2 -(*-l) 



Equations (39) and (40) fulfill the generalized Fourier relation 



J = 



<5(z) dz ' 

where Z) denotes the diffusion coefficient. 

Analogously we define the generalized diffusion velocity v D (z) 



v D (z) 



(39) 

-smh\f(z)-f(a)] 
(40) 

(41) 



1 

w(l-2*) 2 


i 

cosh[/ ( Z ) - /(a)]- (1 - 2*) 2 sinh|> (x) - /(a)] 


_ _ 





rc(z) 



(1 - 2k) 2 cosh[/- (x) - /(a)]- sinh|/(x) - f(a)] 



(42) 

Assuming constant cross section for heat carriers scattering 8(z) = 8 we obtain from 
formula (38) 



/(z) = (l-2*) 2 z, 
/(0)=0, 



(43) 



f(a) = (l-2k) 2 a 
and for density p(z) and current density j(z) 
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j(z) = 



2qwe 



-(l-2k) 2 aS 



1 + Pe-Q-Wos 
(I- 2k) 



(I- 2k) 



-cosh 



(l-2k) 2 -(k-1) 



(2k-l) 2 (x-a)8 



-sinh 



(l-2k) 2 -(k-1) 



(2k-l) 2 (x-a)8 



P(z) = 



2qe 



-(l-2k) 2 aS 



l + /3e 



-(\-2k) 2 a8 



1 



(\-2ky 



-cosh 



(l-2k) 2 -(k-1) 



-sinh 



(2k -l) 2 (x-o) 



(2fc-l) 2 (;t-a)(5 

(l-2fc) 2 -(*-l) ' 
We define Fourier's diffusion velocity Vf(z) and diffusion length, L 



Considering formulae (28) and (29) one obtains 

i 

v F (z) = w(\-2k) 2 , 



L = 



1 



mfp 



1 1 

8(\-2k) 2 (l-2k) 2 
where X mfp denotes the mean free path for heat carriers. 
Considering formulae (44), (45), (46), (47) one obtains 



j(z) = 



2qwe L 



1+ jSe" 



(I- 2k) 



(l-2k) 2 -(k-1) 



-cosh 



(x - a) 



(I -2k) 



sinh 



(l-2k) 2 - (k-1) 



x - a 



P(z) = 



2qe L 



l + /3e~ L 



(i-2ky 

i 

(l-2k) 2 -1 



-cosh 



x-a 



1 



-sinh 



x-a 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(l-2k) 2 -(k-1) 

In Figs. 1-3 we present the results of the calculation for density currents of heat carriers. 
Figs. la-3a are the solution of the Fourier equation for the boundary condition (36). 
Figs. lb-3b represent the solution of the Boltzmann equation formula (48) for a » L,a = L 
and a« L respectively when k e [0,0.5]. For the length of strand, a» L the both solutions 
Fourier and Boltzmann equations overlap. For a < L the Boltzmann equation gives the 
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different description of the transport processes. In that case the solution of the Boltzmann 
equation depends strongly on the scatterings coefficient) of the carriers. Recently [12], the 
heat conduction in one-dimensional system is actively investigated. As was discussed in 
papers [12] the dependence of density current on L can be described by the general formula 

./ ~ - 

where a can be positive or negative. The same conclusion can be drawn from the calculation 
presented in our paper. In this calculation coefficient a depends on the scattering cross 
section for the heat carriers. 

4. Conclusions 

In this paper the Boltzmann equation for heat transport in one-dimensional systems is 
proposed and solved. It is shown that heat current density for a < L strongly depends on the 
mechanism of the heat carrier scattering. The obtained results show that: (a) Boltzmann 
equation allows for the description of heat transport for different ratio alL. (b) The heat 
current density strongly depends on the scattering mechanism for the heat carriers, and on the 
length of the strand. 
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Figure captions 

Fig 1(a) The solution of the Fourier equation for the Cauchy boundary condition 
(formula (28)). (b) The solution of the Boltzmann equation for a » L. 
Fig. 2(a) The same as in Fig. 1(a). (b) The solution of the Boltzmann equation for a = L. 
Fig. 3(a) The same as in Fig. 1(a). (b) The solution of the Boltzmann equation for a « L. 
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1 (a) FOURIER EQUATION 
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2 (a) FOURIER EQUATION 
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3 (a) FOURIER EQUATION 
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